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ABSTRACT 



Satellite-derived ocean surface velocity vectors in the California Current System 
(CCS) are compared with in situ hydrographic and Doppler data. The in situ data 
were acquired during the April 19S1 phase of the Coastal Ocean Dynamics Experiment 
(CODE) experiment. In general, the satellite-derived velocities agreed with the in situ 
data. Due to the baroclinic nature of the study region, the satellite vectors were found 
to be representative of the subsurface geostrophic How. Although the number and 
concentration of the satellite vectors was small, these vectors were capable of resolving 
the mesoscale features located in the study region. Comparison of co-located data 
revealed that the in situ Doppler velocity measurements were approximately 1.5 times 
larger than the satellite-derived velocity vectors. These results agree with a similar 
study by Flament et al. (1985) which was conducted in an offshore region. 
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I. INTRODUCTION 



As the cost of oceanographic hydrographic surveying (i.e. ship time, equipment, 
personnel, etc.) continues to rise, satellite remote sensing of the ocean environment has 
become increasingly popular. Although remote sensing does have some limitations 
(e.g. cloud coverage in Sea Surface Temperature (SST) measurements), satellites do 
allow oceanographic data to be collected over very large areas in an extremely short 
amount of time. This synoptic view of a large region is not routinely available in 
oceanographic research other than by satellite. Satellite sensing also allows researchers 
and ordinary users to "resample" their particular area of interest sometimes several 
times a day. This availability of data has provided, and will continue to provide, the 
impetus for research into determining how to acquire the maximum amount of 
information from the image(s). A particular technique for extracting satellite 
information is addressed in this thesis. Of importance for the Navy is the ability to use 
satellite-derived data in ocean model initialization for prediction purposes. Significant 
advances could also be made in weather forecasting using an air-ocean coupled model 
initialized with such satellite data. 

This research centers on comparing satellite data with in situ hydrographic and 
Doppler data. It is the object of this paper to verify coastal ocean flow vector fields, 
that have been generated using a feature tracking technique on three sequential 
infrared (IR) Advanced Very High Resolution Radiometer (AVHRR) images, in the 
California Current System (CCS). How well these vectors represent the coastal ocean 
flow will be determined by comparing them with in situ data collected during the 
Coastal Ocean Dynamics Experiment (CODE) experiment. Figure 1.1 shows the study 
region. 

Chapter II presents background information on the general character of the CCS 
and previous investigations that have addressed satellite-derived flow vectors. The 
techniques used to generate the velocity vectors and the associated error budget are 
described in Chapter III. Chapter IV compares satellite results obtained in the study 
with hydrographic and Doppler data. Conclusions and recommendations are presented 
in Chapter V. 
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II. BACKGROUND 



A. PHYSICAL OCEANOGRAPHY OF THE CALIFORNIA CURRENT SYSTEM 

Hickey (1979) describes the climatological characteristics of the currents which 
comprise the CCS; the California Current, the California Undercurrent, the Davidson 
Current and the Southern California Countercurrent. The California Current is defined 
as a surface equatorward (low, while the California Undercurrent is defined as a 
subsurface poleward flow over the slope. The Davidson Current is defined as a 
poleward flow that occurs during the fall and winter north of Point Conception. The 
Southern California Countercurrent is also defined as a poleward flow. However, it 
occurs south of Point Conception and inshore of the Channel Islands in the California 
Bight. Since the research in this paper centers on the central California coastal region, 
the current of particular interest is the California Current. 




Fitzure 2.1 Seasonal profiles of the meridional California Current 
oil northern California (from Ikeda and Emery, 19N4). 

The seasonal variation of the California Current is shown schematically in Tigure 
2.1 (Ikeda and Emery. 19S4). There are three axes of strong southward flow. '1 he 
nearshore flow is most fully developed in the spring and early summer south of Cape 
Blanco. The offshore flow on the other hand, is most fully developed in the late 
summer or fall. South of Cape Mendocino, a second offshore region of strong 
southward flow is present throughout the year. The nearshore flow is driven by the 
longshore component of the wind stress, while the two offshore flows are driven by the 
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curl of the wind stress (Hickey, 1979). Various counterllows occur north of Point 
Conception near the coast. During the fall and winter these counterllows displace or 
in some cases even replace the nearshore southward current. While the large scale 
climatology of the California Current can be depicted quite easily, the daily variations 
in its mesoscale conditions are much more complex. 

To investigate these mesoscale features, a field program (CODE) was designed. 
The basic design of the survey grid (Huyer and Kosro, 1987) was influenced by 
previous studies of other upwelling regions and as a result, the alongshore spacing 
between sections was much larger (20-30 km) than the onshore-offshore spacing (5-10 
km). It should be noted that, because of this spacing, dynamic topography could not 
resolve eddies or meanders with alongshore scales less than 40-60 km. Each section 
extended about 45 km from the coast (spanning the continental shelf and upper slope) 
and each of the sections were chosen to be locally perpendicular to the shelf isobaths. 
Conventional CTD casts to a maximum depth of 1000 m and continuous operation of 
a Doppler acoustic log (DAL) to a depth of 150 m (to obtain velocity profiles) were 
included in each section. The geopotential anomaly (AD) of the sea surface relative to 
100 dbar was computed for each survey. For the few stations that were shallower than 
100m, AD 0 ax was calculated directly and ^^p max 'ioo was calculated by linear 
extrapolation from the next offshore pair of stations to obtain AD t) / 100 . 

Huyer and Kosro (1987) found that the resulting geopotential anomaly maps 
were very similar to maps of near-surface salinity. Lowest values of geopotential 
anomaly coincided with the high-salinity upwelled water near shore, and the highest 
values were associated with the low-salinity water offshore. The geostrophic flow 
appeared to be along isohalines, with the low-salinity water on the right and generally 
equatorward over the shelf and upper slope. During the April 1981 survey, a baroclinic 
jet was detected, but its position varied widely over the study region. The strongest 
portion of this jet was located just offshore near Point Arena. 

Because of the continuous operation of the DAL, the current vectors generated 
from the Doppler acoustic log show substantially increased mesoscale variations in the 
flow pattern as compared to the geopotential patterns. Huyer and Kosro (1987) state 
that maps of the current vectors can resolve eddies or meanders with alongshore scales 
of at least 20-30 km (the distance between sections). Although the dynamic 
topography of the April 1981 survey indicates a weak and meandering baroclinic jet, 
the associated DAL maps do not show clear evidence of a continuous equatorward 
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current. Rather, during the April 19S1 survey smaller-scale eddy-like features, which 
penetrated down to at least 150m. were observed over the slope and shelf. A more 
complete discussion of the April 19S1 survey results and a comparison with satellite 
data are presented in Chapter IV. 

Modeling studies (Ikeda and Emery, 1984) of similar mesoscale circulation 
features off the coast of Vancouver Island have shown that the growth of these feature 
is due to the baroclinic instability associated with the vertical shear between the surface 
current and the undercurrent. Ikeda and Emery propose that the mesoscale features 
seen in satellite imagery off the California coast are due to the same mechanism 
because of the great similarity in the mean current structure found off Vancouver 
Island and California. After modifying the nonlinear numerical model for the 
California coast. Ikeda and Emery were able to successfully simulate the initial 
excitation, growth, and cascade of energy to longer wavelengths of these features. 

Flament et al. (1985) studied the evolution of an upwelling filament offshore 
along the coast of California during June 19S1. Velocity vectors were generated by 
tracking the cores of small eddies which Flament et al. assumed to be advected by the 
large-scale flow. A comparison of the satellite-derived velocity vectors against in situ 
DAL current vectors revealed that the Doppler log speeds were larger by a factor of 
approximately 1.5. Flament et al. proposed that either the speed of tracked jet 
increased during the hydrographic survey, or there was a systematic bias of the 
satellite-derived velocities. 

B. REMOTELY SENSED OCEAN DATA 

In the past. AVHRR data has been used to intercompare with in situ sea surface 
temperature measurements (Bernstein, 1982), in situ dynamic height measurements 
(Rienecker et al., 1985), and most recently, to track ice movements (Ninnis et al., 
19S6). Feature tracking research can be broken into several different approaches based 
on determining; 1) how to track features that appear in AVHRR images; and 2) which 
features to track. 

Saunders (1973) used an airborne IR radiometer to estimate the SST of an area 
90 km X 90 km off the New England coast. This area was remapped over a two-week 
period which provided enough data to allow both time derivatives and horizontal 
gradients of SST to be computed. Using this information, Saunders developed a 
technique to compute surface flow. However, Saunders had to make the following 
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assumptions: that the mixed layer temperature was a quasi-conservative quantity, and 
that horizontal diffusion of temperature could be neglected for scales of 5 km and time 
intervals of a few days. The technique predicted that the component of How parallel to 
the mean horizontal temperature gradient was solely due to the rate of change of the 
SST. Unfortunately, the technique could not predict the component of flow normal to 
the temperature gradient unless the horizontal flow in the upper ocean was assumed to 
be horizontally non-divergent. The normal velocity component was then computed by 
determining the normal displacement of an isotherm between a pair of flights. Due to 
a lack of in situ surface current measurements, Saunders was unable to determine 
exactly how accurately this method represented true conditions. Yet, this technique 
represents a start in trying to determine whether surface flow information could be 
obtained from a remote sensing platform. 

Another method was used by Kelly (1983) to compare satellite-derived SST maps 
with in situ current velocity measurements in the CCS. Kelly's method required 
"decimation" of the images by averaging over 5x5 boxes of data. The "decimation" of 
the data consisted of computing EOF's (Empirical Orthogonal Functions) using a 
computationally faster SVD (Singular Value Decomposition) method. However, the 
SVD method had a weakness in its inability to handle missing data and required that 
any missing data (due to cloud coverage) be filled in. These data gaps were filled in 
using a linear combination of data in space and time, which represented the best 
estimate in a least-squares sense. Only "small" gaps were filled by spatial interpolation 
due to Kelly's requirement that a missing datum could only be replaced if the mean 
square error of the estimate was less than 30% of the variance at that location. 
Composite images were created using partial images separated by 12 or 24 hours to fill 
in large holes. The SST EOF's were then computed using these interpolated and 
composite images. Kelly showed that the EOF's of the SST field revealed that the 
shapes of the SST mean and anomalies were dominated by local effects (coastal 
topography, orographic effects, and spatial variations in the wind). The EOF's also 
revealed that the intensity of the anomalies varied with the strength of the local wind. 
Then Kelly showed that the SST patterns were the same as those of the near-surface 
velocities, which suggested that the near-surface current field may have also been 
forced by the local variations in the wind. 

Vastano and Reid (1985) used an interactive algorithm to estimate the motion of 
sub-mesoscale features. This interactive algorithm required the operator to subjectively 
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generate the displacement vectors of individual features in successive pairs of AVI I RR 
images at 12 to 24 hours intervals. The displacement vectors were generated by 
tracking perturbations that had at least a four degree difference with respect to the 
surrounding sea surface temperature. Vastano and Reid proposed that these 
perturbations (which were small relative to the baroclinic radius of deformation) were 
created by small-scale wave or turbulent processes and that their displacement 
represented the local flow regime. Vastano and Reid stated that Rossby-like or 
internal gravity wave motion could be present, but that the wave phase speeds would 
be much less than maximum theoretical values and therefore assumed that these would 
move slower than typical mesoscale surface currents. An additional assumption that 
Vastano and Reid made was that diffusion, heat exchange with the atmosphere and 
vertical motions were much less than the exchange of heat induced by advection over 
the time period analyzed. The generated vectors were then used to compute 
streamfunctions which allowed the nonuniform vector distributions to be transformed 
into a uniform vector field. Vastano and Reid have shown that vectors generated in 
this manner agree quite well with estimated sea surface topography in the area of the 
Oyashio Front northeast of Japan. However, the primary disadvantage of this method 
is that it requires cloud-free views of the sea surface. 

In an effort to move away from the subjectiveness of feature tracking, Emery et 
al. (19S6) developed a method that was objective and operator-independent in 
computing surface advection velocities from sequences of satellite images. In this 
method the flow vectors were generated by identifying the end points of the 
displacement vectors as determined by the maximum cross correlation between a 32 x 
32 pixel search window in the second image and a 22 x 22 pixel template in the first 
image. Once the satellite images were transformed into maps of SST gradient, the 
maximum cross correlation (MCC) technique was applied. The resulting vectors were 
spatially coherent and exhibited excellent agreement with tracks of shallow-drogued 
drifters. However, near coastal boundaries there was a problem in that some of the 
windows contained land and would not allow consistent representations of the 
advective SST changes to be generated. As a result, all such velocity vectors were 
edited out of the final vector fields. Another problem that was detected with the MCC 
method was that it could not respond to rotational or deformational changes while the 
feature-tracking technique could. 



15 



As shown above, each method has demonstrated that satellite-generated velocity 
vectors can be used to determine surface advection velocities, and yet each method also 
has some limitation. To what extent these techniques could provide accurate ocean 
surface velocities in the coastal region is unknown. Therefore, while acknowledging the 
subjectiveness of the feature-tracking technique, this method was chosen because no 
one to date has thoroughly examined what ocean motion these vectors represent in the 
coastal ocean regime. Additionally, the programming required to develop such a 
method was straightforward and easily implemented. 
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III. TECHNIQUE 



For this research AVHRR data (in the form of ocean flow vectors) was 
compared against in situ Doppler current measurements and hydrographic data which 
were acquired during the April 1981 phase of the Coastal Ocean Dynamics Experiment 
(CODE) research off the coast of California. The results of the CODE research (Iluyer 
and Kosro. 1987) were used as ground truth for all subsequent comparisons. The 
satellite data was acquired at Scripps Institution of Oceanography and those passes 
which covered the area of interest during the period of the CODE research were used 
to create the ocean flow field. The April 19S1 time frame was chosen after viewing the 
available satellite images and determining; 1) which images had the least amount of 
cloud contamination in the CODE research region and 2) which images coincided with 
one of the periods examined by H uyer and Kosro. 

A. SATELLITE IMAGERY 
1. Description 

The satellite images that were used to generate the ocean flow vectors were 
512 X 512 pixel portions of three consecutive NOAA-6 orbits. The data were 
calibrated in aperture brightness temperature and resampled onto a so-called "equi- 
rectangular" map grid (which was equivalent to creating a small area plotting sheet 
using a "plane-sailing approximation"). The error associated with this assumption will 
be addressed in a later section. The images were loaded onto the Naval Postgraduate 
School (NPS) COMTAL system and examined to determine the extent and location of 
any cloud coverage which would contaminate the infrared signal the satellite received. 
The result of the comparison between the visible (VIS) and the infrared (1R) images for 
the first pass (04; 27, SI 16:15:00 GMT) was that some thin high clouds were detected 
near the western edge of the CODE region (Figure 3.1). For the middle (night) pass 
(04 2S 81 03:34:00 GMT) some low stratus along the coast was detected by noting the 
changes in the coastline (Figure 3.2). Additionally, some mid and high clouds were 
also present, and these limited some feature tracking. Comparison of the VIS and IR 
images for the third pass (04 28/81 15:52:00 GMT) revealed only a very small amount 
of low stratus along the western edge of the research area (Figure 3.3). 
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Figure 3.1 XOAA-6 AVHRR Infrared image 04 27 81 16:15:00 GMT. 
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.2 NOAA-6 AVHRR Infrared image 04' 28 81 03:34:00 GMT. 
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Figure 3.3 



NOAA-6 AVHRR Infrared image 04 28 SI 



15:52:00 GMT. 
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2. Error Analysis 

The plane-sailing approximation can be described as centering an equi- 
rectangular plane tangent to the earth at a known point in the image(s). For all three 
passes the tangent point was chosen to be Point Reyes, with a known position of 
3S°X, 123°\V. Once the center point of the plane had been selected, the image was 
then projected onto the equi-rectangular grid of the plane. To determine the error 
associated with this approximation, twelve landmarks were selected and their line and 
pixel locations determined. The line and pixel locations were then converted to image 
latitudes and longitudes based on location of Point Reyes in each image. These image 
latitudes and longitudes were then compared against the true latitudes and longitudes 
which were determined using U.S. Department of Commerce Sectional Aeronautical 
Charts for the region of interest. The differences between the image and the true 
latitudes and longitudes revealed two important pieces of information. 

First, all the images needed to be rotated approximately 8 degrees to the left 
(counter-clockwise) in order to align each image with true North. Unfortunately, this 
correction could not be made on the COMTAL due to software limitations (the 
COMTAL would only allow image rotations in 90 degree increments). However, this 
correction was applied through hand calculation and all the coastal landmarks then 
plotted within 2 pixels of their true location. This residual error (approximately 7 - 10 
cm sec) can be traced to the geometry between the satellite, the view angle, and the 
earth location of the pixel(s). While rotating the images eliminated some of the 
navigational error due to applying the plane-sailing approximation, the remaining error 
of 7 - 10 cm'sec cannot be corrected without solving for the associated rotational effect 
on the different spherical triangles involved in the problem. Unfortunately, these 
calculations could not be made because a program, (now being developed at the Naval 
Postgraduate School) which navigates an image by obtaining a best fit for several 
landmarks via changing the satellite's roll, pitch, and yaw, was not available at the time 
this study's AVHRR data were processed. (For more information on this program see 
the recommendations section of Chapter V.) Therefore, it was necessary to reduce the 
image navigation error via the feature tracking program (Oceantrak) as explained in 
the next section. 

The second important piece of information that was expected but could not be 
determined precisely was the East, West error due to using an equi-rectangular grid 
when lines of longitude on the earth converge as the distance from the equator 
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increases. Since Point Reyes was selected as the center point in processing these 
images, the assumed value of 1.1 km per pixel would be valid for any pixel along the 
3S 1 * 1 parallel of latitude. North of this latitude, the "width" of each pixel should get 
smaller by approximately 15 meters per degree of latitude. Similarly, south of 3S°X 
each pixel should get larger by approximately the same amount (15 meters) per degree 
of latitude. However, in the region of this study (Point Reyes to Point Arena) the 
difference in pixel size would only change the computed velocities by a maximum of 
0.037 cm sec. 

B. OCEANTRAK 
1. Description 

A modified cloud tracking program (Oceantrak) was used to generate the flow 
vectors. Oceantrak required the user to input a "crossing" time for each pass so that a 
scaling factor cm ps (centimeters per pixel per second) could be computed. The 
resulting factor was then used to calculate a velocity based on the displacement of a 
submesoscale feature seen in both images. 

The displacement vectors were created by "marking" a feature with a movable 
cursor on one image, and then "toggling" to the second image and "marking" the same 
feature again. Because a water vapor correction was not made to any of the three 
images, the gray shade of most "marked" features changed between images. To avoid 
"marking" different features, unique bends, bumps, and other small scale perturbations 
were chosen to be tracked. Offshore these types of features were quite easily tracked, 
unless cloud cover moved into the area. But, as the coast was approached, it became 
increasingly difficult to locate the same feature from one pass to another, and therefore 
the number of vectors generated decreased as the coast was approached. Additionally, 
features that did not appear to move were tracked by "marking" the center of the 
feature in both images. However, only a few of these types of features were present in 
the images that were examined. 

Oceantrak was also used to reduce the extraneous movement of features and 
landmarks due to an initial lack of image to image coregistration. Because successive 
satellite paths are not the same, the images must be aligned (coregistered) so that 
landmarks will coincide spatially on the screen. Since neither of the sets of passes were 
coregistered, the middle pass was chosen to serve as ground truth for registering the 
other two passes. The images were registered by first determining the movement of 
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several land features and then attempting to minimize the movement in the area of the 
CODE data. By using Oceantrak to coregister the images most of the error from the 
plane-sailing approximation was eliminated. Because the images could only be moved 
by integer pixel amounts, there was some slight residual movement (error) between 
images. 

2. Error Analysis 

For this experiment the satellite's "crossing" time over San Francisco Bay (to 
the nearest minute) was determined. The maximum error resulting from determining 
the "crossing" time to the nearest minute was 0.004 cm p s (less than 0.2 percent). 
While the error associated with the "crossing" time was quite small, the error due to 
imperfect coregistration of the images ranged from 2.7 cm sec from the East for the 
first set of vectors to 3.5 cm sec from the Northeast for the second set. These last two 
errors are due to the limitation created by only being able to move the images by an 
integer pixel amount. 

C. CODE DATA 

1. Description 

The CODE Doppler profile data was acquired through the continuous use of 
an Ametek-Straza Doppler acoustic log (DAL). The horizontal velocity vector profiles 
spanned the 15-150 m depth range and the data points were vertically separated by 6.5 
m. The currents were estimated by filtering the combined DAL and navigation 
(LORAX-C) data over 30 minute intervals. The CODE survey grid was designed to be 
small enough to be completed in 2-3 days, and yet large enough to cover 
approximately 100 km of coastline. Since the time required to complete a survey was 
substantially less than the typical duration (7 to 10 days) of a strong wind-driven 
upwelling event, each of the mesoscale surveys was considered "synoptic", or at least 
quasi-synoptic. 

2. Error Analysis 

For this study the Doppler current vectors were obtained from the Huyer and 
Kosro (1987) paper. The errors associated with these vectors were determined to be 
2-3 cm; sec in the alongshore component and 4-5 cm, sec in the cross-shore component 
(Kosro, 19S5). That portion of the signal due to inertial motion or internal tides is not 
quantified in Huyer and Kosro (1987). The potential influence this submesoscale 
variability might introduce was addressed in another survey (Kosro and Huyer, 1986). 
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After conducting surveys on two jets off the coast of northern California. 
Kosro and Huyer (19S6) concluded that the offshore jets were likely continuations of 
the associated nearshore jets. To determine the submesoscale variability (caused by 
internal tides and inertial motion) in the DAL current and hydrographic measurements. 
Kosro and Huyer separated the (expected) geostrophic mesoscale flow from the smaller 
scale and more rapidly fluctuating components. This was accomplished by using 
objective analysis to estimate the geostrophic current at each DAL vector location 
beyond the 500 meter isobath, and then both sets of vectors were smoothed using a 
Gaussian filter. Once the small scale fluctuations had been removed. Kosro and Huyer 
found that the correspondence between the geostrophic and the DAL currents was 
"striking". Assuming that the above jets were truly extentions of the associated 
nearshore jets, then ignoring internal and inertial motion in the nearshore region 
should not introduce any large errors. For this reason I have not included them in the 
error budget table. 



D. TOTAL ERROR BUDGET 

The largest source of error in this study was the use of an inaccurate navigation 
routine in navigating the images. This method of navigation and the resulting 
imperfect coregistration of the two images introduced errors of up to 2.7 cm sec in the 
computation of all velocity vectors generated using the first and second images. 
However, the resulting imperfect coregistration from using this method of navigation 
introduced errors ranging as large as 3.5 cm/sec in the generation of the second set of 
velocity vectors. While these errors are large (in comparison to the error in the pixel 
size due to the converging/diverging line of longitude), a qualitative comparison 
between the Doppler data and the satellite-derived vectors is still be possible since the 
error in the DAL data is as large or larger. See Table 1 for a breakdown of the 
different error components. 
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TABLE 1 

TOTAL ERROR BUDGET 



ERROR 

SOURCE 


NORTH - SOUTH 
ERROR (cm/sec) 


EAST 

ERROR 


Plane Sailing 
Approximation 


. 00023-. 00025 


0. 034 


Crossing Time 


0. 07 - 0. 08 


0. 07 


Imoerfect 
Coregistration 
Pass 1-2 


0. 00 




Im.Derf ect 
Coregi st rat ion 
Pass 2-3 


2. 48 




Doppler 
Acoustic Log 
( ?Cosro , 196o ) 


2.0 - 3.0 


4. 0 



- WEST 

( cm/sec 

- 0. 037 

- 0. 08 

. 70 

. 48 
-5.0 



IV. RESULTS 



As expected from California Current climatology, the general flow in the CODE 
region was towards the south during the April 19S1 survey. The dynamic topography 
as determined from the data acquired during the April survey also agreed (in general 
terms) with the climatological mean dynamic topography for the month of April 
(Hickey. 1979). The winds during the April survey were strongly favorable for 
upwelling (Huyer and Kosro. 19S7). 

A. CODE SURVEY 

All of the following descriptions of hydrographic and DAL features are my 
interpretations of data presented in Huyer and Kosro (19S7). 

1. Hydrographic Features 




Figure 4.1 Maps of in situ sea surface temperature and temperature 
at 20 m (from Huyer and Kosro. 1987). 

Plots of the April survey SS T showed the coldest surface water adjacent to the 
coast, and a strong offshore temperature gradient, both of which are consistent with 
the strong upwelling favorable winds (Figure 4.1). Maps of the surface salinity also 
reinforce the fact that upwelling was occurring during the study period, in that the 
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more saline water was located adjacent to the coast (Figure 4.2). Both the temperature 
and salinity distributions at 20 meters also verify that upwelling was taking place 
during the April survey. The map of geopotential anomaly of the sea surface relative 
to 100 dbar was very similar to the map of the near-surface salinity, with the highest 
values being associated with the low-salinity water offshore and the lowest values 
coinciding with the high-salinity water near the coast (Figure 4.3). Huyer and Kosro 
propose that southward advection of subarctic water would explain the persistently low 
surface salinities offshore, while northward advection of local runoff through a narrow 
and sporadic countercurrent along the coast would explain the isolated pockets of low 
salinity water over the shelf. Huyer and Kosro also state that the circulation pattern 
they observed during the survey was much more complex than they had expected. In 
fact, on the northern edge of the survey area, all of the above maps indicate the 
presence of a warm, low-salinity, anticyclonically rotating eddy. A vertical phase shift 
of this eddy can be seen in the temperature maps (see Figure 4.1). The surface position 
of the eddy as shown by the SST map agrees well with the satellite imagery. The 20 
meter position of the eddy (which is southwest of the surface position) as shown by the 
temperature plot at that level agrees with the plot of the DAL data at the 20 meter 
level. In the southwest region of the study area there is little indication (due to the 
spacing of the hydrographic stations) of a second warm, low-salinity, anticyclonically 
rotating eddy present in the above maps. Unfortunately, a comparison of the surface 
and 20 meter positions of this eddy cannot be made due to the coarse resolution of the 
hydrographic data in the area of interest. 

2. DAL Features 

The 20 meter map of the DAL current vectors (Figure 4.4) for the April 
survey shows the southern eddy quite well, but does not give as good a representation 
of the eddy to the north. The DAL map at 20 meters also shows a general southward 
flow in the region. 

At 80 and 150 meters, the DAL current velocity maps (Figure 4.4) are very’ 
similar to the 20 meter map. However, as expected, the number of near shore velocity 
vectors decreases as the water shoals. Offshore, both of the eddies are represented to 
the same degree as they are in the 20 meter map and the general southward flow of the 
region is also still evident. The only real differences between the 20, 80. and 150 meter 
maps are (as expected) that the magnitude of all the velocities decreases as the depth 
increases. 
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Fieure 4.2 Maps of in situ surface salinitv and salinity at 20 m 
(from Huyer and Kosro, 1'9S7). 




Figure 4.3 Map of the eeopotential anomaly (m 2 s 2 ) 
of the sea 'surface relative to 1 10~dbar (from Huyer and Kosro, 19S7). 
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Fieure 4.4 Maps of DAL current vectors at depths of 20. SO. and 
150 m (Irom Huyer and Kosro. 19S7). 

B. SATELLITE DATA 

In the study region, the one feature that stands out quite clearly is the large 
warm eddy along the northern edge (figure 4.5). Although this eddy was partially 
covered by clouds in the middle (or night) pass, the majority of the feature was 
detectable, particularly the boundaries. Unfortunately, the second eddy to the south 
was not as detectable. In fact the southern edge of this feature was quite diffuse and 
as a result, a boundary could not be determined from the SST imagery. This 
homogeneity of the coastal water mass was also a problem in generating near-shore 
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velocity vectors. The near-shore problem was compounded by coastal fog that 
developed during the night and severely limited use of the middle pass in generating 
velocity vectors near the coast. 

Outside the study area (off Monterey Bay) there were several distinct SST 
features that were easily tracked and which appeared in all three images used in this 
study. Velocity vectors were generated in this area in the hope that another in situ 
data set could be found to serve as a ground truth. Unfortunately, 1 could not find 
such a data set. These other vectors were left on the images so that the reader could 
see how fast some of the features moved, and more importantly how active the coastal 
region is on a mesoscalc basis (Figures 4.5 and 4.6). 

The SST features in the study region, as shown by the second set of satellite 
images, became quite diffuse. This is due in part to the fog which developed during the 
night and obscured many of the nearshore features. Another factor is that a large 
cloud mass can be seen moving into the study region from the northwest. This cloud 
mass and its associated water vapor also contaminated a portion of the satellite sensed 
data for the region. Because of these factors, approximately half as many velocity 
vectors were generated using the second set of images as were generated using the first 
set (Figure 4.6). 

C. COMPARISON OF CODE AND SATELLITE VECTORS 

Because only ten vectors were generated from the second set of images, the 
following comparison will concentrate on those vectors generated from the first set of 
satellite images, where nineteen vectors were generated. Figure 4.7 shows the first set 
of satellite vectors (large arrowheads) superimposed on the 20 meter DAL 
measurements. 

1. Agreement 

Overall, the first set of satellite velocity vectors that were generated in the 
study area agree quite well (in direction of How) with the DAL vectors (Figure 4.7). 
Where comparisons could be made (both types of data located in approximately the 
same position and the directions of flow agreed), the DAL measurements were 
approximately 1.5 times larger than those of the associated satellite derived velocity 
vectors. Eight out of eleven vectors agreed with DAL measurements where the two 
data sources coincided. An additional eight satellite vectors were generated at 
locations in the study area, but outside the domain of the DAL measurements so that 
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Figure 4.5 Satellite vectors generated from passes 1 and 2. 
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Figure 4.6 Satellite vectors generated from passes 2 and 3. 
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Figure 4.7 Map of 20 m DAL and satellite-derived velocitv vectors. 

Larue arrowheads are satellite generated vectors, ' 
small arrowheads are DAL measurements. 

comparisons could not be made. The anticyclonic rotation of the northern eddy is 
verified by the satellite vectors, as well as the general southward flow of the region. 
The satellite vectors also show some support for the southern eddy and show rotation 
in the correct direction (anticyclonically). 

Although Huyer and Kosro (19S7) proposed a sporadic countercurrent along 
the coast could explain the isolated pockets of low salinity water over the shelf, the 
DAL measurements do not reveal any countercurrent. Similarly, the satellite-derived 



vectors Jo not indicate the presence of any countercurrent along the coast. In fact, 
nearshore both the DAL measurements and the satellite generated velocity vectors 
agree quite well. 

2. Disagreement 

There are three satellite-derived velocity vectors which disagree with the DAL 
measurements (Figure 4.7). The first is due to frontal waves propagating southward 
along the sharp boundary between the cold upwelled water and the warmer offshore 
water. These waves exist along the northern edge of the study region (where they are 
ISO degrees out of phase with the northward DAL data). One of the resulting velocity 
vectors (3 cm sec, which is near the noise level of the technique) is present in Figure 
4.7 This vector is located northwest of Pt. Arena, over the 1000 meter isobath. The 
frontal waves also follow the western boundary of the northern eddy (where two 
additional waves were tracked). When these features were originally tracked, they were 
so clear in the imagery that the eddy appeared to be rotating cyclonically. It was only 
after determining that the features were not part of the eddy but rather traveling along 
the western boundary that separated the warm offshore water from the cold upwelled 
water that 1 realized they had to be frontal waves. 

The second disagreement between the satellite vectors and the DAL 
measurements occurred along the northern edge of the southern eddy. In attempting 
to track the diffuse edge of the eddy, one northward vector was generated which was 
approximately 180 degrees out of phase with the southward DAL vector. This vector 
is located southwest of Pt. Arena, again over the 1000 meter isobath and is labelled 20 
cm s. This northward shift in the location of the eddy is most likely due to 
atmospheric contamination from a small group of clouds that moved into area during 
the night. 

There is a third vector that disagrees with the DAL measurements, but the 
velocities of both vectors are small (4-5 cm/s). This satellite vector is located over the 
150 meter isobath, practically due west of Pt. Arena and is labelled 4 cm/s. The vector 
was generated by tracking what appeared to be a perturbation on the edge of the warm 
eddy. While it is possible that the perturbation was an instability on the eddy, clouds 
had moved into this particular area between the two passes used to generate the vector. 
It is possible that the perturbation is atmospheric in nature, not oceanic. 
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3. Other Physical Processes 

Because the winds during the April 1981 survey were upwelling favorable, one 
might expect significant offshore directed velocity vectors associated with Ekrnan 
transport. An estimate of this contribution to the satellite derived motion can be 
obtained from simple analytic solutions (see Csanady, 1982). Csanady's solutions 
provide cross-shelf velocities which are zero at the coast and increase exponentially 
offshore. The maximum cross-shore response occurs at a distance offshore equal to 
the external Rossby radius of deformation (350 - 500 kilometers). However, Kelly 
(1983) has shown that coastal topography plays a significant role in the wind forced 
response of the circulation in the CODE region. The wind driven ocean response is 
thus spatially variable. The satellite vector fields (Figures 4.5 and 4.6) show no 
obvious uniform offshore velocity trend. Based on Csanady's analytic solutions, strong 
offshore velocities are suspected to occur further offshore. 

The good agreement between the satellite-derived vectors and the subsurface 
geostrophic flow indicates that significant contributions by inertial motion are not 
present in the satellite vectors. This is also in agreement with the results of Kosro and 
Huyer (1986), as discussed above. 
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V. CONCLUSIONS AND RECOMMENDATIONS 



A. CONCLUSIONS 

Satellite-derived ocean surface velocities have been compared with in situ 
observations. In general, the satellite-derived velocities agree with the hydrographic 
surveys and the Doppler acoustic log data. This implies that the vectors during this 
time period are representative of the subsurface geostrophic flow’. Additionally, the 
satellite vectors do not show any signs of any inertial motion or Ekman forcing. 

The satellite-derived velocity vectors agreed with the hydrographic data that was 
acquired during the April 19S1 CODE survey. But, while the in-situ hydrographic data 
could not resolve the mesoscale jet and eddies located in the study region, the satellite 
generated vectors did resolve these features. The satellite vectors were found to be 
representative of the subsurface geostrophic flow, due to the baroclinic nature of the 
study region. While this technique can generate velocity vectors which are 
representative of the geostrophic flow, it is important to note that the success of this 
technique depends on the water vapor content of the atmosphere. In particular, the 
effects of fog (prevalent during during upwelling conditions) may preclude the use of 
this technique. 

Like the offshore study by Flament et al. (1985), DAL measurements were 
approximately 1.5 times larger than the satellite-derived velocity vectors (when they 
were co-located). However, the number and concentration of the satellite vectors were 
substantially smaller than the DAL data. This sparseness of data made point by point 
comparisons of the DAL versus satellite data difficult. 

B. POTENTIAL USES 

There are many potential uses for this technique, both as currently demonstrated 
and with some minor improvements. For example, if the concentration of velocity 
vectors could be increased, via either feature tracking from color images or by applying 
an atmospheric correction (so that isotherm displacements might be tracked), the 
technique could definitely provide hypothetical input data for global ocean circulation 
modeling and possibly enough data to initialize a regional ocean circulation model. In 
its current operating condition, the technique could be of assistance in SAR (Search 
And Rescue) operations by allowing the user to visualize ocean surface flows (direction 
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and speed) in any area of interest. This information could then he used to organize the 
search to take maximum advantage of assets (aircraft and or ships) available in 
locating the victim in distress. Information on current variability even - 6-12 hours 
would also be quite valuable in tracking pollution (oil spills) and would assist in 
predicting if. when and possibly where, such pollution might come ashore. The lead 
time this information would provide could then be used to notify local residents and 
position clean-up equipment in those areas where it would do the most good. This 
technique could also be used in ASW (Anti-Submarine Warfare) by allowing features 
of interest (fronts and eddies) to be tracked and their location and movement briefed to 
those forces in the associated region. On a more scientific note, this same information 
would be very useful in assisting shipboard hydrographic surveys of both fronts and 
eddies by providing the ship(s) with close to real-time (shortly after the satellite 
image(s) have been processed) location and movement of such features and by allowing 
gaps in the in situ data to be filled via the information obtained from the velocity 
vectors. The above are just a few of the potential uses of this and similar techniques: 
undoubtly more uses and more accurate information will be found as research in this 
field continues. 

C. RECOMMENDATIONS 

In order to obtain more accurate information from this technique, the satellite 
images that are used must be navigated to the highest degree possible. A new AVHRR 
navigation program is presently being developed at the Naval Postgraduate School (Dr. 
James L. Mueller personal communication, 1987). This program minimizes the error in 
the location of several user selected landmarks by adjusting the satellite's roll, pitch, 
and yaw. Once the overall error has been minimized, the image is resampled and 
remapped with a latitude and longitude computed for every sixteenth line and pixel. 
This approach to navigating satellite images, together with image resampling onto a 
conformal map projection (i.e. transverse Mercator) would significantly reduce the 
distortion which was the dominant source of error in determining displacement vectors 
in this study. Considering the above, it is recommended that if additional research is 
conducted on this technique that the AVHRR navigation program be used when it 
becomes available. 

An additional recommendation for further research with this technique is that 
color enhanced images might provide a denser concentration of velocity vectors by 
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possibly allowing smaller temperature differences to be tracked. However, to color 
enhance the satellite images properly, an atmospheric correction will have to be made 
to allow an accurate determination of the SST. This correction would be necessary to 
avoid color shifts between images that are only the result of an increase or decrease in 
water vapor in the atmosphere. Since this correction requires two 1R channels, only 
satellite data from NOAA-7 and odd numbered follow-on satellites can be used. 
Fortunately, there are additional surveys (during other periods) contained in the Iluyer 
and Kosro (19S7) report to allow a comparison to be conducted using color tracking 
(NOAA-7 images) and in-situ data. 

There are two potentially fruitful areas of study which only require one IR 
channel's worth of data. The first recommendation is to investigate a no wind or low 
wind situation to determine the extent of influence the wind has or does not have on 
the concentration and magnitude of the satellite-derived velocity vectors. Again, the 
Huyer and Kosro (1987) report contains information acquired during three surveys in 
which the winds were quite weak. This would allow additional satellite generated 
velocity vectors to be compared with in situ hydrographic data. The second 
recommendation is to investigate an offshore region in an effort to determine if 
satellite-derived velocity vectors agree or disagree with in situ hydrographic data. 
Although a study similar to the above recommendation has been mentioned (Flament 
et al., 19S5), a detailed comparison between in situ data and satellite generated velocity 
vectors in a coastal transition region has not been published as far as this author 
knows. 

Finally, a word of advice, to maintain an oceanographic frame of mind, avoid 
generating vectors just for the sake of generating vectors and determine if each vector 
is consistent with the local thermal contrast in SST and geostrophy. By maintaining 
the above thought process, the user will avoid generating false vectors associated with 
large scale wave propagation. It should be noted that the above reasoning will provide 
good guidance in generating velocity vectors in most cases. When salinity becomes as 
important as temperature in determining the density, the above reasoning would no 
longer apply. This type of situation, in which salinity appeared to be at least as 
important as temperature in determining surface density, was encountered by Kosro 
and Huyer (1986) in an offshore study. Therefore, any potential users of this or other 
vector generating system should have a good working knowledge about the 
oceanographic conditions that should be expected in the region where they are 
applying this technique. 
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